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ABSTRACT: 


The  design  formulae  and  theoretical  performance  of  a 
maximum-likelihood  estimator  for  towed-array  flow-speed  estimation 
are  discussed.  Under  appropriate  conditions,  it  is  shown  that  a 
maximum-likelihood  estimate  of  convection  velocity  (which  is  related 
to  tow  speed)  is  achieved  by  performing  a cross  correlation 
between  the  outputs  of  two  streamwise-separated  flush-mounted 
surface  probe  hydrophones.  The  estimate,  then,  is  derived  by 
seeking  the  time  delay  which  maximizes  the  cross  correlation. 

Conditions  are  examined  consistent  with  achieving  accurate 
speed  determination  to  within  a fraction  of  a knot.  Fortuitously, 
the  method  derived  has  tlie  property  that,  as  speed  is  increased, 
both  the  level  and  bandwidth  of  the  pressure  signals  increase, 
thus  providing  increased  signal-to-noise  ratio  and  more  accurate  time- 
delay  information. 

It  is  suggested  that  a confidence-level  measurement, 
comparing  the  space-time  correlation  function  derived  from  the  data 
with  a theoretical  model,  could  be  used  to  indicate  flow  anomalies. 
These  flow  anomalies  might  be  attributable  to  disturbances  such  as 
flow  separation.  In  such  case,  the  device  could  be  used  also  to 
sense  flow  stability. 

Two  implementations  are  discussed:  one  using  a DElay  Line 
Time  Compressor  (CELTIC)  and  serial  processing;  the  other  using 
Shift  Registers  (SR)  and  parallel  processing.  A trade-off  of 
serial  f aster-than-real-time  CELTIC  processing  speed  requirements 
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I 
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against  the  high  degree  of  parallel  processing  needed  with  the 

SR  approach  is  discussed.  | 

A reconmiendation  is  made  to  design,  build,  and  test  such  | 

a speed  estimation/stability  sensor  on  standard  towed  arrays.  It 
further  is  suggested  that  empirical  results  and  experience  thus  | 

derived  could  be  extrapolated,  along  with  a study  of  available 
other  array  surface  pressure  coherency  data,  to  determine  the  * 

feasibility  of  building  such  a sensor  for  non-standard  towed  array  | 

application. 


1.0  INTRODUCTION: 


The  report  discusses  the  design  and  performance  of  a 
maximum-likelihood  estimator  for  towed-array  flow-speed  estima- 
tion. It  is  shown  that,  under  appropriate  conditions,  the 
maximum-likelihood  estimator  consists  of  a correlator  which 
measures  the  time  delay  corresponding  to  the  peak  of  the  broad- 
band correlation  measured  between  the  outputs  of  two  streamwise- 
separated  flush-mounted  surface  probe  hydrophones. 

The  hydrophones  have  a circular  sensing  surface,  a 
small  fraction  of  an  inch  in  diameter,  and  respond  primarily  to 
the  high-wavemomber , convective  pressures  produced  by  the  turbu- 
lent boundary  layer.  The  probes  are  separated  longitudinally 
by  a fraction  of  an  inch,  so  that  over  most  of  the  operating  band 
the  coherence  is  high. 

This  method  of  flow-speed  estimation  is  capable  of 
accurate  speed  determination  to  within  a fraction  of  a knot.  It 
has  the  desirable  feature  that,  as  speed  is  increased,  both  the 
level  and  bandwidth  of  the  pressure  signals  increase,  providing  in- 
creased signal-to-noise  ratio  and  more  accurate  time-delay  information. 

If  disturbances  such  as  flow  separation  exist  in  the 
vicinity  of  the  surface  probes,  it  is  expected  that  erratic  results 
will  be  obtained.  Therefore,  it  is  possible  that  this  device  also 
can  be  used  to  indicate  flow  stability.  However,  this  subject 
requires  further  study. 
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It  should  be  pointed  out  that,  prior  to  use  at  sea,  the 
speed-measurement  system  must  be  calibrated  in  a water  tunnel. 

The  calibration  is  expected  to  depend  to  some  extent  upon  the 
location  of  the  sensors  within  the  towed  array  (since  macroscale 
properties  of  the  flow  may  change) . In  addition,  since  the 
relation  between  convective  and  flow  speeds  will  change  with 
speed  (for  a fixed  system  operating  bandwidth),  it  will  be  neces- 
sary to  calibrate  the  system  over  the  entire  operating  speed 
range . 

Originally  this  system  was  conceived  for  application  to 
standaid  towed  arrays.  However,  it  may  be  possible  that  a 
similar  system  can  be  designed  for  speed  and  flow-stability 
estimation  on  other  arrays.  The  technology  exists  today  for 
building  such  sensors.  It  is  li)cely  that,  with  some  types  of 
towed  arrays,  the  turbulent  boundary  layer  will  exhibit  greater 
sensitivity  to  tow-point  and  flow  dynamics  than  a standard  array. 

In  addition,  the  report  derives  the  estimation  accuracy 
to  be  theoretically  anticipated  from  a maximum- liJceli hood  estimate 
of  convection  velocity.  The  errors  perturbing  the  estimation 
accuracy  are  treated  in  two  categories:  a)  bias  errors;  and  b)  fluctu- 
ation errors.  The  bias  errors  are  residual  errors  causing  a 
departure  of  the  average  of  the  estimate  of  the  convection  velocity 
from  its  true  value.  The  fluctuation  error  is  one  of  zero  expecta- 
tion whose  rms  value  measures  the  dispersion  of  the  estimates 
around  their  mean.  It  is,  in  part  due  to  background  noise 
corrupting  the  estimate.  Additionally,  in  as  much  as  the  Turbulent 


Boundary  Layer  (TBL)  pressure  fluctuations  are  noiselike,  a 
residual  fluctuation  error  occurs  when  the  estimate  is  performed 
over  a finite  number  of  independent  samples. 

The  report  considers  the  functional  dependence  of  these 
errors  on  the  space-time  correlative  nature  of  both  the  TBL  and 
background  noise  fluctuations,  signal-to-noise  ratio,  and  the 
number  of  independent  samples  used  in  estimating  the  correlation 
function.  Design  formulae  indicating  the  dependence  of  the 
precision  on  each  of  the  variables  are  developed.  In  addition, 
timing  and  time- resolution  requirements  are  established  consistent 
with  the  number  of  independent  samples  needed  and  the  required 
precision . 


Two  fundamentally  different  implementations  are  discussed: 
one  using  a DElay  Line  Time  Compressor  (CELTIC)  and  serial 
processing;  the  other  using  Shift-Registers  (SR)  and  parallel 
processing.  It  is  shown  that  a clocking  rate  of  the  order  of 
hundreds  of  MHz  is  needed  with  the  CELTIC  to  cover  the  whole 
speed  range.  The  SR  implementation  leads  to  a far  more  modest 
clocking  rate  at  the  expense  a high  degree  of  parallel  processing 
(of  the  order  of  100  channels) . However,  it  is  suggested  that  the 
SR  implementation  is  compatible  with  state-of-the-art  LSI  technology. 

The  report  also  discusses  a scheme  for  deriving  a 
measure  of  confidence  level  to  be  attributed  to  the  estimate  of 
convection  velocity.  This  conf idence -level  measure  is  derived  by 
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comparing  estimates  obtained  from  measured  data  with  stored 
models  of  the  TBL  correlation  structure.  As  such,  a low  con- 


fidence level  would  be  indicative  of  a radical  departure  from 

the  models.  This  departure  then  would  indicate  the  presence  of 

flow  anomalies.  Moreover,  comparison  of  confidence  levels  derived  ^ 

with  a non-standard  array  against  similar  measurements  with  a 

standard  array  could  indicate  the  adequacy  of  extrapolating  the  f 

TBL  models  to  the  non-standard  array  geometry;  wherein,  greater 

sensitivity  to  tow-point  and  flow  dynamics  might  be  anticipated.  ^ 

The  report  recommends  that  available  surface-pressure  | 

coherence  data  for  non-standard  arrays  be  studied  to  determine  the 
feasibility  of  building  a speed-estimation/flow-stability  sensor  I 

for  application  to  them.  Some  of  the  experience  necessary  for  this  r 

step  could  be  acquired  through  measurements  made  with  a similar  ^ 

sensor  applied  to  standard  arrays.  j 
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2.0  CONCLUSIONS  AND  RECOMMENDATIONS: 


It  is  concluded  that  an  estimate  of  the  TBL  convection 
velocity  may  be  made  with  suitable  precision  by  performing  a 
Maximum  Likelihood  Estimate  (MLE) . The  MLE  requires  that  a 
sampled-data  estimate  of  the  cross  correlation  between  the 
forward-  and  the  aft-hydrophone  TBL-pressure  signals  be  formed. 
In  forming  this  cross  correlation,  the  time  series  derived  from 
the  forward  hydrophone  is  time-delayed  to  bring  it  to  the  point 
of  maximum  correlation  with  the  aft  hydrophone. 

In  order  to  improve  the  stability  of  the  estimate  of 
convection  velocity,  the  TBL  pressure  fluctuations  being 
a random  process,  multiple  samples  are  needed.  For  the  rms 
value  of  the  fluctuations  around  the  average  estimator  value 
to  improve  by  being  inversely  proportional  to  the  square- 
root  of  the  nur.iber  of  samples,  N,  these  time  samples  should  be 
selected  far  enough  apart  in  time  so  as  to  be  decorrelated.  Two 
partially  correlated  time  series  are  being  cross-correlated. 
Therefore,  the  requirement  to  achieve  the  full  l/'^N  improvement 
is  that  the  space-time  correlation  coefficient  of  the  TBL  be 
small  for  the  hydrophone  separation,^,  selected  and  for  a total 
time  separation  equaling  the  time  between  samples  minus  that 
time  delay, T,  which  maximizes  the  cross  correlation.  A 5-milli- 
second interval  between  scunples  is  suggested  for  free-stream 
velocities  down  to  8 kts. 


It  further  is  shown  that  the  rn>s  error  of  the  estimator 


i 


is  inversely  proportional  to  an  equivalent  angular-frequency  band- 
width, [(2nf^)^  + (1/3)  (ttW)  ^ , and  the  hydrophone  separation  , C . 

Here,  W is  the  observation  bandwidth,  centered  around  a frequency 
fg.  In  fact,  this  improvement  is  offset  by  another  factor  equal  to 
unity  plus  the  exponential  of  1.4  times  the  Strouhal  number  Cf^/U^, 
E,  being  the  fore-and-aft  separation  of  the  hydrophones  and  the 
convection  velocity.  As  a result,  a minimum  fluctuation  about  the 


average  of  the  estimator  is  shown  to  occur  at  a separation. 


^opt ' 


equal  to  1.584  hydrodynamic  wavelengths  (i.e.,  to  1.584  U^/Iq) . This 
restriction  enforces  separations  of  the  order  of  5x10  ^m(~0.2") 


at  free-stream  velocities  of  the  order  of  16  kts. 


Additionally  it  is  shown  that  the  estimator  rms  fluctua- 
tions decrease  with  the  Signal-to-Noise  (Power)  Ratio  (SNR),  y,  being  ' 

directly  proportional  to  (l+y  *)•  Moreover,  the  average  of  the 
estimator  exhibits  a slight  bias  error  (i.e.,  its  statistical  average  ^ 

is  not  equal  to  the  true  value)  for  low  SNR  and  a small  number  of 
samples.  Typically,  an  SNR  of  100  is  adequate  to  remove  the  bias 
error  and  ensure  a negligible  contribution  to  the  flucturation  error. 

However,  for  the  fluctuation  error  to  be  acceptably  small,  say  1% 
of  the  maximum  velocity,  at  least  1000  independent  samples  would  be 
needed.  These  would  ta)<e  5 seconds  to  accumulate  and  this  determines 
the  time  constant  of  the  velocity  sensor  in  response  to  changes  in 
velocity.  The  1000  independent  samples  are  shown  to  be  more  than  r 

adequate  to  make  the  bias  error  insignificant. 

» 
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It  IS  apparent,  from  the  theoretical  analysis,  that  the 
MT.F  of  convection  velocity  provides  an  estimation  accuracy  exceed- 
ing that  implied  by  the  Cramer-Rao  bound.  This  anomaly  appears 
to  be  real,  and  probably  has  its  roots  in  the  non-regular  nature 
of  the  TBL  space-time  function  assumed;  which  is  attributed  to 
Bakewell  [2].  Indeed,  Cramer  [3]  shows  an  example  where  non- 
regularity can  lead  to  abnormally  high  precision.  Swerling  [4] 
has  pointed  out  similar  anomalies,  and  associates  them  with  a 
tigliter  bound  acquired  through  a procedure  due  to  Barankin . (cf 
Swerling  loc.  cit.  and  Glove  [5] ) Be  that  as  it  may,  the  MLE 
does  lead  to  a useful  and  intuitively  satisfying  procedure  for 
estimating  the  convection  velocity.  The  errors  predicted  also 
have  a heuristic  explanation.  Therefore,  it  is  suspected  that  it 
is  the  Cramer-Rao  bound  which  is  not  tight  enough  for  this  non- 
regular  case,  and  not  that  there  are  any  deficiencies  with  the  MLE 
procedure . 

It  is  recommended  that  in  designing  a velocity  sensor, 
only  a two-to-one  range  of  velocity  be  spanned  for  one  fixed  hydro- 
phone-pair  separation  and  that  the  optimum  separation  of  1.584  U^/f^ 
be  chosen  for  the  geometric  mean , ■^^rdnu\  max  U^.  , of  the  velocity 
range . 

To  limit  errors  to  values  of  the  order  of  1%  up  to  16  kts , 
the  time-delay  increments  should  be  of  the  order  of  10  micro- 
seconds. Two  implementations  are  discussed:  one,  using  a faster- 
than-real- time  serial  processing  scheme  centered  around  DElay 
Line  Time  Compressors  (DELTIC) ; the  other,  a real-time  parallel 
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scheme  using  Shift  Registers  (SR).  The  DELTIC  scheme,  consistent 
with  the  10-microsecond  real-time  resolution,  necessitates 
clocking  rates  of  the  order  of  several  hundred  megahertz.  For 
this  reason,  this  approach  is  not  recommended.  By  way  of 
contrast,  the  SR  scheme  requires  a 200-stage  SR,  with  the  last 
100  stages  tapped,  running  at  a shift  rate  of  100  KHz.  Beyond 
that,  100  parallel  1-bit  multipliers  and  100  parallel  10-bit  up/ 
down  counters  (reset  every  5 seconds)  would  be  needed;  again 
operating  at  a clocking  rate  of  100  KHz.  Modern  LSI  technology 
easily  can  accommodate  this. 

The  SR  approach  can  be  made  compatible  with  an  LED,  4- 
digit,  floating-point  display  of  free-stream  velocity.  In  addi- 
tion, it  can  be  made  to  provide  a feature  for  displaying  a 
confidence  level.  This  would  facilitate  providing  the  display 
operator  with  an  indication  of  any  anomalous  behavior  of  the 
axial  flow  along  the  towed  array;  for  example,  a velocity  sensor 
near  the  tail  of  the  array  might  be  a good  indicator  of  flow  and, 
hence,  array  stability. 

To  derive  a confidence- level  measure,  models  of  the 
clipped  space-time  correlation  appropriate  to  each  of  the 
convection  velocities  associated  with  each  correlator  time-delay 
tap  could  be  stored  in  a Read  Only  Memory  (ROM) . The  model 
appropriate  to  the  estimated  velocity  then  would  be  accessed  and 
compared  through  differencing  either  in  a mean-square  or 
absolute-value  sense  summed  across  the  time-delay  taps  dispersed 
either  side  of  the  tap  for  the  estimated  velocity.  A small 
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value  of  this  measure  would  correspond  to  a high  confidence 
level,  say  9 on  a scale  of  0 to  9.  Lower  numbers  corresponding 
to  higher  values  of  the  measure  would  be  an  indication  that  the 
estimated  cross-correlation  function  did  not  match  well  with  the 
mode 1 . 


It  is  recommended  that  a velocity  sensor/flow-stability 
sensor  be  designed,  built  and  tested  on  standard  towed  arrays. 

The  resultant  data  could  be  reduced  and  compared  with  the  theory 
over  a wide  range  of  tow  speeds  and  sensor  locations.  It  further 
IS  recommended  that  available  surface-pressure  coherence  data  for 
non-standard  arrays  be  studied  to  determine  the  feasibility  of 
building  a speed-estimation/flow-stability  sensor  for  applications 
to  other  than  standard  arrays. 
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TECHNICAL  DISCUSSIONS: 


In  attempting  to  estimate  the  convection  velocity  associated 
with  flow  in  the  Turbulent  Boundary  Layer  (TBL) , it  is  assumed  that 
measurements  will  be  made  with  a pair  of  f lush“mounted  hydrophones 
separated  along  the  longitudinal  flow  direction  of  a towed 
array.  Each  hydrophone  should  be  small  in  its  longitudinal 
spatial  extent  in  order  that  sufficient  time  resolution  of  the 
pressure  fluctuation  might  be  achieved;  see  Corcos  [1].  Moreover, 
their  spatial  separation  should  be  small  enough  to  preserve  some 
degree  of  space-time  correlation.  In  this  respect,  assuming  that 
broad-band  pressure  fluctuation  measurements  will  be  made  consistent 
with  good  time  resolution,  a model  attributed  to  Bakewell  [2]  will 
be  presumed. 

One  type  of  hydrophone,  which  would  respond  primarily  to 
high-wavenumber,  convective  pressures  produced  by  the  TBL,  could  have 
a circular  sensing  surface.  In  order  to  sense  flow  distortion  when 
an  array  angle-of-attack  is  invo’.ed,  three  sets  of  hydrophone  pairs 
spaced  at  0°,  120°,  and  240°,  respectively,  around  the  circumference 
of  the  array  appear  to  be  appropriate.  For  speed  sensing,  each 
hydrophone  pair  should  be  located  in  a region  of  (quasi)  spatially 
homogeneous  and  temporally  stationary  TBL  flow  over  the  towed 
array.  In  a location  near  the  tail-end  of  the  array,  anomalous 
results  caused  by  flow  instabilities  might  be  anticipated.  Thus,  the 
velocity  sensor  also  might  be  used  as  a detector  of  flow  instability. 

It  appears  appropriate  to  consider  performing  a Maximum- 
Likelihood  Estimate  (MLE)  of  the  convection  velocity  in  lieu  of 
attempting  a Baysian  Estimate  (BE) . The  BE  requires  a priori  and 
(error)  cost  information  not  readily  available;  whereas  the  MLE  can 
be  estimated  assuming  only  knowledge  of  the  appropriate  multi- 
dimensional joint  Probability  Density  Functions  (PDF's).  The  MLE 
maximizes  the  probability  of  a correct  decision  (estimate),  without 
regard  for  an  incorrect  decision  (estimate)  or  its  cost.  However, 
in  the  absence  of  a priori  and  (error)  cost  information,  the  MLE 
represents  the  most  acceptable  alternative.  In  order  that  a suffi- 
ciently smooth  estimate  of  the  convection  velocity  may  be  derived, 
multiple  time  scimples  are  needed;  gathered  at  time  intervals  sepa- 
rated enough  to  provide  uncorrelated  samples  in  a tim.e  series 
extracted  from  each  hydrophone  of  the  pair.  Of  course,  correlation 


3.1 


■r 


must  be  maintained  between  the  time  series  if  the  convection 
velocity  is  to  be  deduced.  In  this  study,  consistent  with  low 
implementation  complexity  and  expense,  spatial  smoothing  achieved 
by  using  more  than  two  hydrophones  will  not  be  considered. 

In  order  to  investigate  the  form  of  the  MLE , it  will  be 
assumed  that  the  small  flush-mounted  hydrophones  measure  the 
pressure  fluctuations  dominated  by  the  convective  energy  of  the 
TBL.  The  only  corrupting  teirni  to  be  considered  will  be  broad-band 
extraneous  noise.  Inasmuch  as  the  noise  may  be  caused  by  sources 
which  exhibit  spatial  correlation  over  the  small  distance  separating 
the  hydrophone  pair,  this  possibility  will  be  considered  in  deriving 
the  MLE. 


The  broadband  nature  of  the  signals  caused  by  the  TBL 
pressure  fluctuations  also  the  noise  fluctuations  favors  using  a 
real  rather  than  an  analytic  signal  representation.  The  observed 
time  series  derived  from  the  forward  hydrophone  in  the  TBL  flow  stream 
is  given  by 
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+ N 


likewise,  from  the  aft  hydrophone 
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,t 


The  quantities  x 
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x<^),t 
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and  X 
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x(^>,t 


are  the  spatial  locations  along  the 
axial  direction  of  the  center  of  the  forward  hydrophone  and  aft 
hydrophone  respectively.  The  quantities  t=t^^^  and  t=t are  the 
times  at  which  the  time  samples  are  derived  from  each  of  the  forward 
and  aft  hydrophones  respectively.  Indeed,  N time  samples  are  derived 
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This  notation  will  be  contracted  to 


in  each  of  the  time  series  at  times  t 
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for  the 


for  the  aft  hydrophone 


x(f’)  = p(F)  ^ 
and 


x(A)  = p ; 

where  the  bar  under  each  quantity  indicates  that  it  is  a column  vector 
N-tuple,  and  a superscript  T transposes  a column  vector  into  row 
vector.  Then, 
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Under  the  assyu^tion  that  both  the  TBL  pressure  and  noise 
fluctuations  are  Gaussianly  distributed  stochastic  processes,  the 
above  formalism  leads  to  a 2N-variable  joint  PDF  conditional  on 
the  convection  velocity  given  by 


det 


(F ) 

In  the  PDF  equation,  the  column  vector  X is  adjoined  as  a 

(A) 

subvector  to  the  other  sub-vector  X so  that 


j,(A) 


also  the  bar  over  X indicates  a statistical  ensemble  average 
{expected  value)  which  is  assumed  to  be  zero  for  real  representations 
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of  the  broad-band  stochastic  processes.  The  matrix  R is  the 
correlation  matrix  which,  because  of  the  vanishing  averages  of  X, 

IS  the  same  as  the  covariance  matrix.  The  matrix  R ^ is  the 
inv('rse  of  R also  det  (R)  is  the  determinant  of  R.  The  correlation 
matrix  R is  given  by  the  expected  value 


This  matrix  has  symmetric  form.  The  submatrices  in  the  R matrix 
are  given  by  the  diagonal  matrices 


0 


0 


where  the  diagonal  form  results  from  choosing  the  time  samples 
sufficiently  far  apart  so  as  to  decorrelate  them  on  both  an  intra- 
and  interspatial-channel  basis.  Only  the  cross-correlation  matrices 
T T 

x(F)  and  X^^*  X^F)  ^^e  functions  of  the  convection  velocity, U^. 


Invoking  the  previously  mentioned  spatial  homogeniety  and 


temporal  stationarity , and  exploiting 
submatrices,  leads  to 

the  diagonal 

forms  of 

the 

w (x|u  ' 

. (x„ 

y - 20x>'n 

zN  \ ' C ^ 

where 

' n=l 

- 

r = = 

= '^x' 

9 

r = = 

= °x' 

also 

Y (f’)  ' 

Xi 

. (A)  _ 

^1  “ 

V (F)  y (A) 

■ »x  ®x’ 

X 

X/F)  = 

= X,/F) 

X.  = p o ^ • 

^ ^X  °x 

1 

1 

Again  it  should 
velocity,  U^. 

be  noted 

that  only  is  a 

function  of  the  convection 

I 

I 

I 
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4 


should  be 


At  this  point,  the  general  form  of  and 

evaluated.  Under  the  spatial-homogeniety , temporal-stationarity 
assumptions,  for  any  n=l, ,N, 


= (p 
\ n 


(F)  ^ „ (fA 

n / 


likewise 


Moreover , 


Px^x^  = X X 

XX  n n 


.(p  >'■>  * N <^>)  (p  <''>  + N <^>) 
\n  n/\n  n/ 


= PpOp^  + 


Therefore , 


Pv  = 


‘•p°p’  ^ Vn’ 


“p“  * 
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In  order  to  derive  an  MLE  of  the  convection  velocity,  U^,  a 
solution  is  sought  to  the  equation 


-9  log^ 


where  log^  is  the  natural  logarithm.  This  equation  leads  to 


I (x  + [x  - 2p  X ^^^X 

9U_  ,\n/  \n/  ^Xn  n 

c n=l-‘» 7 r 


+ N log 


(^-Px' 


Performing  the  partial  differentiation  and  factoring  out  leads  to 

) ° c 

This  leads  to  two  solutions  appearing  as  the  second  and  third  central 
equations;  either 


N n=l  L \ I -J  ' ' 


The  second  central  equation  is  not  data-dependent . There- 
fore, it  only  indicates  that,  given  a priori  the  form  of  the  space- 
time  correlation,  the  correct  procedure  to  derive  an  MLE  of  is  to 

find  the  value  of  U where  the  correlation  reaches  a maximum.  A 

c 

minimum  is  not  appropriate  for  the  broad-band  case.  The  third  central 
equation  is  data-dependent. 


The  term  i I (x 

Nn=l'^"  / 


is  a consistent  estimator  of  (i.e.. 


as  N-*®  it  converges  in  probability  to  However,  according  to 

Cramer  [3],  it  is  not  an  unbiased  estimator  of  This  biased 

_ ^ I fx  ... 


1 Wx  ‘f’)’ 

property  may  be  corrected  by  replacing  the  term  \ n / by 

1 ^ / (F)\^ 

'(N=TT  ^ V^n  / ' ^ consistent  and  unbiased  estimator  of 

n=  1 ' 

, (F ) 

, assuming  the  expected  value  of  X is  zero.  Likewise,  by 


introducing  a multiplication  factor 


the  estimators  of  the 


variance  of  and  the  covariance  of  X^^^X^^^  may  be  rendered  both 

consistent  and  unbiased.  Thus,  the  third  central  equation  may  be 
replaced  by^^ — 

A 2 ^ 

Further,  it  is  reasonable  to  introduce  a„u,  and  p„i„  as  the  consistent 

X I N X I N 

and  unbiased  estimates  of  the  variance  and  covariance,  namely 


xIn  (n 


-1)  rn  ) (N-1)  ) 


n n - ^ r X X . 

Px|N°X|n  ■ (N-1)  n n 
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(F)  (A) 

again  under  the  assumption  that  the  expected  value  of  X and  X 
are  zero.  In  this  way  the  third  central  equation  may  be  compacted 
in  the  form 


^ 2 

^ ^ 2 2 (N- 1 ) I - 

^XlN^XlN 

• 4 

The  expected  value  of  the  third  central  equation  is 

-N  A,  2 

2 {N-l)Oj^^- 

^x1n°x|n 

(N-1) 

f ^ 1 \ ^X  X ' 

r ^ 

using  the  fact  that 

\ / 

^X|N°X|N^ 

Px^x'- 

The  only  way  in  which  this  equation  can  be  satisfied  exactly 
IS  if 


that  is  if 

Px  = -1 


or 


in  the  case  of  broadband  fluctuations  only  the  first  condition  is 
appropriate.  When  this  condition  applies  the  fourth  central  equation 
becomes 

A ” ' 

^X|N®XjN  ^ °X^ 

This  condition  implies  that 


which  is  consistent  with  the  first  central  equation. 

The  condition  for  an  MLE  of  then  is  to  apply  a (usually 
gradient)  search  loop  to  locate  the  maximum  of  the  estimate  Pj^jjg  of 
the  TBL  space-time  correlation.  In  a bandwidth  W arithmetically 


-4 


I 


centered  around  a frequency,  f q , Bakewell  [2]  has  assumed  that  the 
TBL  space-time  correlation  is  given  by  a fifth  central  equation 
(shown  in  Figure  1) 


(,  (C,t)  = exp|-0.7  -^n 

sin  ttw| 

1 

^ L “c  J 

X cos  ^2tT  f Q - 

■n  W 

[-(5/u,); 

1 

here  the  Strouhal  number  is  given  by 


I (A)  (F)  1 

where  5 is  the  axial  separation,  |x  -x  ),  of  the  hydrophone 
centers,  and  i is  the  time  delay  between  the  corresponding  time  Scim- 
ples  in  each  of  the  time  series,  i.e,. 


, = t - t = 

T 


^ ^ (A)  _ t (F) 

'"n  '"n 


^ . (A)  (F) 

N N 


At  any  time  delay  satisfying 


" = U~ 
c 


the  TBL  space-time  correlation  achieves  a local  maximum 


= exp  -0. 7 \Qf 


= exp[-0.7|T|fJ 


Obviously,  if  the  condition  for  the  third  central  equation  is  to  apply 
exactly,  the  space-time  correlation  of  the  TBL  pressure  fluctuations, 
for  a spatial  separation  C and  a time  delay  t = C/U^,  must  approach 

closely  to  unity.  This  implies  that  C must  be  very  much  smaller  than 
the  hydrodynamic  wavelength 


C<<X,  = c <<  C 
n — j-  ; 
0 ^0 
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5 m/s 


By  way  of  an  example,  the  hydrodynamic  wavelength  when  = 0.6Uo„ 
(i.e.  , Uoo  =8.3  m/s  = 16.2  kts)  , at  = 1,000  Hz,  is  given  by 
= 5xl0"  ^m  = 5 mm. 

Fortunately,  the  second  central  equation  is  more  tolerant 
with  respect  to  allowing  much  wider  hydrophone  pair  separation.  For 
example,  when  ^ is  equal  to  the  hydrodynamic  wavelength,  Aj^,  also 
assuming  that  the  SNR  and  the  number  of  samples  are  very  large, 
then  pjj  = 0.5.  Under  these  conditions  the  second  central  equation 

is  very  close  to  being  exactly  satisfied.  In  fact,  for  an  SNR  of 
100  and  ^ = 100  tne  left  side  of  the  equation  is  within  1% 
of  the  right  side.  This  arises  from  the  first  central 
equation  yielding  in  the  worst  case,  when  the  noise  is  correlated 
(i.e.  , = 1 ) , 

Pjj  = 0. 99  Pp  + 0.0099 

for  a SNP  of  100  (i.e.,  a power  ratio  of  100).  When  the  hydro- 

phones are  separated  by  a hydrodynamic  wavelength  pp  = 0.5,  corre- 
sponding to  Pj^  = 0.  5049.  Furthermore,  when  N = 100,  for 
p^  = 0.5049  the  third  central  equation  leads  to 

"xllo;  °X  |1PP  . r i ^]p 

= •*  ■ f 0.5049 

[99  99*1.2551 

= 0.502/ 
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■* 


0.5  to  within  0.4%. 


this  approximates  pp  = 


For  the  case  where  the  noise  background  has  a partial 

correlation  coefficient  p.,,  in  order  to  achieve  on  the  average 

N 

relatively  unbiased  estimates  with  error  smaller  than  ±e,  the  SNR 


Y 


also,  the  number  of  independent  samples  needed  for  each  hydrophone 
channel  is 

333. 3 


N> 


l+^l+2eyp2 


- 165.7  . 


A couple  of  examples: 


a.  For  the  worst  case  of  completely  correlated 

noise,  and  a TBL  space-time  maximum  correlation 
of  0.5,  to  achieve  an  average  error  of  1%  or 
less  requires  100  independent  samples  and  a SNR 
of  100.  For  an  uncorrelated  background,  the  SNR 
requirement  remains  at  100.  The  0.5  peak 
correlation  would  be  achieved  with  a hydrophone 
pair  separation  of  one  hydrodynamic  wavelength. 


b.  Again  for  completely  correlated  noise,  and  a 
TBL  space-time  maximum  correlation  of  0.1,  to 
achieve  an  average  error  of  1%  requires  164 
independent  samples  and  a SNR  of  890.  For 
an  uncorrelated  background,  the  SNR  requirement 
reduces  to  100.  The  0.1  peak  correlation  would 
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be  achieved  with  a hydrophone  pair  separation  of 
3.29  times  the  hydrodynamic  wavelength. 


At  this  point,  the  variance  of  the  estimator  should  be  derived 
in  order  to  determine  the  fluctuation  error  around  the  mean  of  the 
estimator.  If  the  estimator  were  a regular  estimation  case  of  the 
continuous  type  (see  Cramer  [3])  the  estimation  variance  could  be 
compared  with  the  Cramer-Rao  bound  to  determine  the  efficiency  of 
the  estimator.  Using  the  definition  in  Cramer  [3] , it  appears  that 
the  estimator  is  not  sufficient;  i.e.,  W2n(— ^^c)  be  factored 

as  the  product  of  ^2N^— '32N^— ^ indepen- 

dent of  the  U^.  The  Cramer-Rao  bound  can  be  satisfied  only  when  this 
condition,  and  the  condition  that 


3 

3 


lo9e 


U 

c 


k 


) 


(where  k is  independent  of  U^  but  may  depend  on  U^)  both  are 
satisfied.  It  therefore  appears  that  the  convection  velocity 
estimator  will  not  be  an  efficient  one,  i.e.,  it  will  not  achieve  the 
Cramer-Rao  bound.  However,  because  of  the  discontinuity  in  Bakewell's 
space-time  correlation  function  at  ^-►O  or  the  regularity 

assumption  is  open  to  suspicion.  Indeed,  one  manifestation  of  non- 
regularity can  be  a lower  variance  of  the  estimation  error  than 
given  by  the  Cramer-Rao  bound  (cf.,  Cramer  loc.  cit,  page  485).  In 
this  case,  another  tighter  greatest  lower  bound  for  all  relatively 
unbiased  estimates,  over  all  cases  (regular  and  non-regular),  and 
all  SNR's  is  given  by  Barankin's  approach;  see  Swerling  [4], 

Glove  [5] . 


In  assessing  the  variance  of  the  error  between  the  estimate, 

A 

U^,  of  the  convection  velocity,  and  its  true  value,  U^,  the  following 

IS  appropriate.  A high  enough  SNR  is  assumed  such  that,  in  the 

vicinity  of  the  peak  of  the  correlation  function,  no  ambiguous 

secondary  peaks  occur.  Moreover,  for  the  purpose  of  estimating  the 

variance,  the  number  of  samples  used  is  assumed  to  be  large  enough 

to  render  the  estimator  relatively  unbiased.  As  previously  indicated, 

the  estimated  value  of  the  convection  velocity  is  given  by  the 

solution  U = U to  the  data  dependent  form  of  the  second  central 
c 

equation,  viz: 

^^X|N  = 0 
9U 

Around  the  value  U = U^,  this  derivative  has  the  mean 

slope 

\3U’  jv  - 

This  is  the  second  (i.e.,  linear)  term  in  a Taylor's  series  expansion 
of  the  first  derivative  of  made  around  Hence,  the  error 

A 

U^-U^  may  be  written  with  good  approximation  as 


The  construction  is  shown  graphically  in  Figure  2.  It  follovs  that 
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the  variance  of  the  error  may  be  written  as 


However,  by  the  theorem  for  differentiating  products 


The  term  on  the  left  vanishes  because  of  the  skew  symmetric  form  of 
the  product  which  averages  to  zero;  also 


9U 


0 


Therefore , 


- -'XIN 


9^0 


X N 


9U- 


Remembering  that  the  correlation  estimator  is  a function  of  the 

given  fixed  separation , C / and  the  time-delay  t,  which  in  turn  is 


equal  to  C/U , it  may  be  written  as  0 


X |N 


(U) 


In  this  way; 


9^lS 


- rx  N- 


X N 


9U' 


U = U. 


\ 9U^  > 


u = u. 


maybe  written  as 
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where  the  covariance  function  is  given  by 


Hence,  the  variance  of  the  error  in  estimating  the  convection 
velocity  may  be  written  approximately  as 


dimensionally,  this  will  be 
^ and  covar  ^^e 


the  square  of  a velocity  as 
numbers . 


Returning  to  the  definition  of  namely, 


(u)  x^<'''(ui, 

n=  1 


in  order  to  highlight  the  functional  dependency  on  £;  as  well  as  U, 
it  may  be  written  in  the  two-dimensional  form 


^J^x(x,t„|u^)  x(xH.VT|U^), 


where  the 
ve locity 

^|N 


conditionality 

is  given  by  U . 

c 

be  expressed  as 


on  U indicates  that  the  actual  convection 
c 

Using  this  notation,  the  covariance  of 
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covar 


It  is  assumed  that,  under  a steady-flow  situation,  the  pressure 
fluctuations  of  the  TBL  sensed  by  the  hydrophones  are  represented 
by  temporally  stationary,  spatially  homogeneous  Gaussian  processes. 
Moreover,  the  broad-band  noise  fluctuations  corrupting  the  observa- 
tions also  are  considered  as  represented  in  a similar  way.  The 
space-time  correlation  (normalized)  structure  of  the  TBL  pressure 
fluctuating  field  is  considered  to  be  in  accord  with  the  fifth 
central  equation.  Furthermore,  the  hydrophone  spacing,^,  will  be 
comparable  with  a hydrodynamic  wavelength,  which,  in  the  range  of 
convection  velocities  under  consideration,  is  very  much  smaller  than 
an  acoustic  wavelength.  I’or  this  reason,  the  corrupting  noise  field 
may  be  considered  as  completely  spatially  correlated.  In  which  case, 
the  normalized  space-time  correlation  coefficient  for  the  noise 
fluctuations  in  a band  of  width  W centered  around  a frequency  f^  is 
given  by 


sin  [ttWt  ] 
I ttWt  ) 


cos  I2irf qT]  . 


Under  the  above  conditions,  it  is  shown  in  appendix  A that 
the  variance  of  the  estimation  error  is  given  by  the  horrendous  but 
dimensionally  correct  and  intuitively  satisfying  equation 
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In  the  case  where  the  SNR  is  very  high  (i.e.,  y ^ becomes  very  small), 
this  central  equation  reduces  to 


In  the  above  central  equation, 


^2TifQ)  " + (1/3)(ttW) 


2 


is  a measure  of 


second  moment  of  the  power  spectral  density  of  the  pressure  fluctua- 
tions. As  such,  it  is  a measure  of  (equivalent  angular-frequency  bandwidth)^. 
It  includes  the  fact  that  a bandpass  region  of  width  W centered  around  f^  is 


considered.  If  f^  is  made  too  large  relative  to  W (i.e.,  narrow  band 
measurements  of  the  TBL  pressure  fluctuations  are  extracted) , the 
space-time  correlation  function  will  exhibit  a discernible  quasi- 
periodic  structure  possessing  ambiguous  peaks  which  nearly  compete  with 
the  central  peak  in  determining  the  extraction  of  an  estimate  of  U^. 

To  avoid  this  situation,  the  bandwidth  ratio,  ^q/W,  of  the  broad- 


band  TBL  pressure  fluctuations  should  be  no  greater  than  3/4. 

This  extreme  situation  illustrating  the  secondary  peaks  of  the 
space-time  correlation  coefficient  is  shown  in  Figure  3.  In  fact,  inas- 
much as  the  lower- frequency  region  of  a towed  array  hydrophone 
output  usually  is  dominated  by  low-wavenumber  energy  induced  by 
vibration  coupling  and  module  scattering,  it  is  prudent,  anyway, 
to  provide  a low-pass  cutoff  at  the  bottom  of  the  band  of  interest. 
Likewise,  beyond  some  upper  frequency  determined  by  the  boundary 
layer  thickness,  a knee  occurs  in  the  power  spectral  density  of  the 
TBL  pressure  fluctuations.  Beyond  this  the  TBL  power  spectral 
density  begins  to  roll  off  at  a rapid  rate.  Then,  above 
this  knee,  even  with  equalization,  the  contribution  of  the  TBL  power 
would  eventually  be  masked  by  the  ambient  noise  contributions. 

Therefore,  it  would  be  imprudent  to  operate  beyond  some  high-pass 
cutoff  frequency  determined  by  the  anticipated  boundary-layer  thick- 
ness. These  low-  and  high-pass  cutoff  characteristics  can  be  made 
consistent  with  the  bandwidth  ratio  requirement  previously  cited. 

The  central  equation  appropriate  when  the  TBL  pressure 
fluctuations  dominate  over  the  ambient  noise  (i.e.,  in  the  high-SNR  case) 
then  indicates  that  the  inaccuracy  due  to  estimation  (fluctuation) 
noise  decreases  as  the  equivalent  bandwidth  increases.  This  may  be 
achieved  by  increasing  both  f^  and  W consistent  with  the  previously 
cited  restrictions.  In  addition,  the  accuracy  improves  with  greater 
hydrophone  spacing  but  is  offset  through  multiplication  by  the  sum 
of  unity  and  an  exponential  function.  The  optimum  spacing  is  given 
by  differentiating  the  quotient  of  the  sum  of  unity  and  the  exponential 
function  in  the  numerator  and  in  the  denominator  with  respect  to  ^ 
and  equating  to  zero.  The  optimum  value  then  is  given  by  a solution 
to  the  transcendental  equation 


This  equation  has  a solution 


C = 1.584U^ 


i.e.,  the  separation  should  be  approximately  1.584  hydrodyneimic 
wavelengths.  In  this  case, 


which,  for  = 3/4,  becomes 

w“ 


0.0896  U 2 

varfU  -U  1 = T5 

\ c cj  N 


or 


0. 3U 


This  implies  that,  to  achieve  a 1%  error  in  the  fluctuations, 
at  least  896  independent  saunples  are  needed  consistent  with  a hydro- 
phone separation  of  ^ = 1.584U  /f 

c 0 


Applying  the  optimum  separation,  of  1.584  hydrodynamic 
wavelengths  and  a bandwidth  ratio,  ^q/W,  of  3/4  to  the  full  equation 
for  tne  variance  of^U^-U^^,  involving  the  noise  as  well  as  the  TBL 
pressure  fluctuation  terms,  leads  to  a close  approximation  for  yi  10 
given  by 
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By  way  of  an  example;  with  1.584  and  fQ/W  - 3/4 


(e.g. , fp  = 1,500  Hz,  W = 2,000  Hz),  to  achieve  a Root-Mean-Square 

(rms)  fluctuation  error  of  0.94%  of  the  actual  convection  velocity 
requires  about  1,000  independent  samples.  With  the  parameters 
chosen,  pp  = 0.33,  and  = 0.05.  From  previous  considerations, 

for  the  bias  error  to  be  less  than  0.34%, 


Ni 


333.3 


165.7 


or 

134 

Obviously,  N = 1,000  more  than  satisfies  this  lower 
bound.  Moreover,  the  required  SNR  consistent  with  a bias  error  of 
0. 34%  is 


y = 100  0.05  - 0.33 

0.34  0.33 

= 250 

In  other  words,  when  the  number  of  independent  samples,  N,  and 
the  SNR,  y,  are  chosen  large  enough  consistent  with  a small 
fluctuation  error,  a negligible  bias  error  is  encountered.  The 
rms  addition  of  the  bias  error  and  the  fluctuation  error  leads 
to  a total  error  of  less  than  ±1%  of  the  convection  velocity. 
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4 . 0 IMPLEMENTATION  CONSIDERATIONS  : 

For  a free-stream  velocity,  , of  16.2  kts,  the  convection 
velocity,  U^,  is  about  5 m/s.  With  set  at  1,500  Hz,  the  hydro- 

- 3 

phone  separation  would  have  to  be  5=5.28  x 10  m = 5.28  mm  (i.e., 

5 = 0.208").  The  band  of  frequencies  being  observed  would  be  500 

Hz  to  2,500  Hz.  In  addition,  the  travel  time  between  hydrophones  in 

- 3 

the  TBL  would  be  about  10  s (i.e.,  1 millisecond).  To  measure  to 
within  better  than  1%  would  require  a time  resolution  of  about  10  ^s 
(i.e.,  10  microseconds).  To  perform  the  processing  digitally  then 
requires  a sampling  rate  of  about  100  KHz.  However,  in  order  to 
decorrelate  the  samples  in  the  time  series  derived  from  a hydrophone, 
they  should  be  selected  less  frequently  than  about  2500  Hz  (i.e.,  at  an 
interval  of  400  microseconds),  say  2000  Hz.  This  would  ensure  that 
intrachannel  samples  did  not  correlate  with  each  other.  However, 
when  the  time  series  extracted  from  the  fo3n«/ard  hydrophone  is  time- 
shifted  relative  to  the  time  series  extracted  from  the  aft  hydro- 
phone (see  Figure  4)  care  must  be  ta)cen  to  avoid  inducing  inter- 
channel correlation  between  other  samples.  To  avoid  this  artificially 
induced  correlation,  the  time  between  samples,  (t  -t  .),  minus  the 

maximum  anticipated  time  delay,  x=^,must  be  sufficiently  large  so 

U 

c 

as  to  make  the  value  of  the  space-time  correlation  function, 

.-j.  ('  ,t^-t^-»-4  I > negligibly  small  for  all  m/n.  For  free-stream 

velocities,  U^,  down  to  8 kts  an  interscunple  period  of  5xl0~^s 
would  suffice.  In  that  case,  in  order  to  accumulate  1,000  inde- 
pendent samples,  an  observation  period  of  at  least  5 s would  be  needed. 
Moreover,  if  a Shift  Register  (SR)  were  used  to  achieve  the  time 
delay  necessary  to  bring  the  forward  hydrophone  into  time  registra- 
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FIGURE  4:  FORWARD  & AFT  TBL  FLUCTUATION  SAMPLE  CROSS  CORRELATION. 


FIGURE  5a: 


EXAMPLE  OF  DELTIC  TIMING  FOR  100  DIFFERENT  TIME  DELAYS. 
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FIGURE  5c:  FORWARD  HYDROPHONE  SIGNAL  DELTIC 


tion  with  the  aft  hydrophone,  at  least  200  stages  would  be  needed 
for  a speed  range  down  to  8 kts  consistent  with  1%  velocity  error 
up  to  16  kts,  A 2-to-l  velocity  range  appears  to  be 
reasonable  without  reoptimizing  the  spacing  between  hydrophones. 

The  spacing  should  be  optimized  for  , "^/ni  n (U^) . max ( U^) . 

Inasmuch  as  both  the  TBL  pressure  and  the  noise  fluctuations  are 
broad-band,  a clipper-correlator  would  not  significantly  degrade 
the  extraction  process.  Either  a serial  CELTIC  correlator  or  a 
parallel  SR  correlator  could  be  used. 

Figure  5 shows  the  correspondence  needed  between  time  samples 
for  the  uLLTIC  iniplementation . Unfortunately,  in  order  to  perform 
all  the  correlation  functions  with  a time-delay  resolution  of 

- 5 

10  S in  real  time,  a time  speedup  of  1000:1  forces  clocking  rates 
in  the  hundreds  of  MHz  region.  Although  specialized  correlators 
have  been  tested  in  the  laboratory  using,  for  example,  a combination 
of  Gun  and  charged  storage  diodes,  in  an  operational  device  their 
use  is  not  recommended. 

The  existing  LSI  technology  would  allow  the  fabrication  of  a 
miniturized  processor  using  multistaged,  tapped  SR's.  A schematic 
for  performing  the  multiple  correlations  shown  in  Figure  4 by  paral- 
lel processing  is  shown  in  Figure  6.  The  number  of  1-bit  multipliers 
and  10-bit  up/down  counters  needed  must  match  the  number  of  tapped 
stages  on  the  SR.  Typically  for  a velocity  range  spanning  a 2-to-l 
ratio,  and  consistent  with  1%  velocity  uncertainty  at  the  highest 
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FIGURE  6:  PARALLEL  CHANNEL  SHIFT  REGISTER  IMPLEMENTATION. 


velocity, an  SR  with  200  stages  and  taps  on  the  last  100  stages 
would  be  required.  Thus,  this  form  of  implementation  does  not 
appeal  to  challenge  today's  LSI  technology. 

In  Figure  6,  logic  for  peak-selecting  (i.e.,  determining  and 
indexing  the  greatest  of  the  100  up/down  counter  states)  is  shown. 

Each  tap  corresponds  to  a different  free-stream  velocity,  so  that 
an  LED  display  could  be  made  to  indicate  velocity  through  a decoder. 

In  addition,  a confidence  check  could  be  provided  by  measuring 
the  mean  squared  deviation  between  the  output  of  the  correlator  on  each 
tap , and  a Read-Only  Memory  (ROM)  stored  set  of  correlation  coefficients 
(i.e.,  a model  of  the  TBL  space-time  correlation  function)  appropriate 
for  the  convection  velocity  corresponding  to  each  tap.  Namely  , the 
ROM  would  store  ^arcsin  [pp(5  T ~ rA) ] corresponding  to  the  clipped 

space-time  correlation  function  centered  at  a particular  time  delay 
f = The  index, r,  would  ideally  span  either  side  of  the 

correlation  peak.  This  would  require  additional  taps  and  stages  on 
the  SR  in  order  to  accommodate  the  whole  2-to-l  velocity  range. 

The  confidence  number  would  be  derived  by  differencing  the 
data-dependent  correlator  outputs  and  the  appropriate  TBL  space-time 
correlation  function  model  correctly  indexed  and  selected  from  the 
ROM  by  the  peak-selector.  These  differences  could  either  be 
squared,  or  their  magnitude  extracted,  then  added  over  the  range  of 
r.  The  resultant  value  would  indicate  the  degree  of  match  between 
the  data  and  the  model  and,  hence,  a measure  of  the  quality  of  the 
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velocity  estimate.  This  confidence  level  could  be  indicated  by 
one  digit  on  the  LED  display. 

What  has  been  described  is  one  possible  implementation  of  a 
convection  (hence,  free-stream)  velocity  indicator.  The  theoretical 
aspects  of  its  accuracy  have  been  analyzed,  showing  feasibility. 
Moreover,  implementation  is  believed  to  be  well  within  the  state 
of  the  art. 
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APPENDICES 


A. 1 Deviation  of  the  Variance  of  the  Estimator  Error 


It  was  shown  in  § 3.0  that  the  variance  of  the  estimator 
error  was  given  by  , 


,3U' 


covar 


var 


fu  -u  'I  = 
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U=U 


3^ 


i3U' 


U=U 


In  order  to  evaluate  the  numerator  and  denominator  of  this  expression, 
the  relationship  for  is  recouched  as 


N 


Z X(x,t^)  X(x+^,t^+x) , 


in  order  to  highlight  the  spatial  and  temporal  dependency.  In  this 
case 


Thus,  in  evaluating  the  numerator  of  the  variance  of  the  estimator, 
the  covariance  may  be  written  as 


-Al- 


Since  all  of  the  X{  , ) functions  comprise  the  sum  of  two  stochastic 
processes,  P(  , ) and  N(  , ),  which  are  both  Gaussianly  distributed, 
then  X(  , ) is  Gaussianly  distributed.  In  such  case  the  average  of 
the  fourfold  product  may  be  replaced  by  the  sum  of  three  terms,  each 
term  being  the  product  of  the  average  of  product  pairs  extracted  from 
the  fourfold  product  in  cyclic  order.  One  of  these  three  terms  cancels 
the  negative  term  appearing  under  the  double  summation  shown  in  the 
last  equation.  It  follows  that 


Remembering  that 


X (x  , t)  X (x+f, , t+T  ) = [P  (x,t)  + N (x  ,t)  1 [P  (x+5  , t+T  ) + N (x+C  , t+T  ) ] 


Op^Pp(C,T)  + 


as  both  the  TBL  pressure  and  the  noise  fluctuations  have  average  zero, 
are  assumed  to  be  temporally  stationary  and  spatially  homogeneous,  and 
are  mutually  uncorrelated,  then  the  covariance  may  be  expanded  as 
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Moreover,  for  5 small  compared  with  an  acoustic  wavelength  (In 
§ 3.0  it  is  shown  that  an  optimum  choice  for  ^ is  1.584  hydrodynamic 
wavelengths;  a hydrodynamic  wavelength,  , being  very  much  smaller 
than  an  acoustic  wavelength,  A , for  low  Mach  number  flow.)  the  space- 
time  structure  of  band-limited  noise  (of  bandwidth  W centered  around 
frequency  f^)  is  given  with  good  approximation  by 


Pfj  (CfX) 
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Here,  C is  the  velocity  of  sound,  and  the  minimum  acoustic  wavelength, 
min  A^,  corresponds  to  the  highest  frequency  in  the  noise  band.  These 

statements  relative  to  the  space-time  structure  of  the  noise  may  be 
interpreted  as  stipulating  that  the  noise  is  highly  spatially 
correlated  for  small  hydrophone  separations , irrespective  of  the 
spatial  distribution  of  the  noise  field;  e.g. , spherically  isotropic, 
cylindrically  isotropic,  etc. — in  fact,  any  noise  not  arising  from  a 
mechanism  similar  to  that  generated  by  TBL  flow.  Generally,  the  noise 
corrupting  the  observations  arises  from  quite  different  mechanisms,  such 
as  surface-wave  noise  or  distant  shipping.  In  such  case,  the  previous 
assumption  about  the  space-time  structure  of  the  noise  is  valid. 
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From  the  above,  it  follows  that 
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Collecting  terms  and  inserting  into  the  numerator  of  the 
covariance  expression  leads  to  a form 
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and 


Returning  to  the  denominator  of  the  expression  for  the 
variance  of  the  error, 


and 
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This  is  the  same  generic  form  as  the  covariance  term.  Therefore, 
the  general  expressions  to  be  doubly  differentiated  wi th-respect- to 


U are 


fj^(U)  = 


sin 

- ^ _ 


~ u ) 

' c 'J 
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These  forms  may  be  reduced  further  to  one  canonical  form  by  substitutions 
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respectively.  The  canonical  form  then  becomes 
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where  u can  take  on  the  value  u^^  or  U2' 
In  either  case. 
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It  then  follows  that 
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This  second  limit  is  a measure  of  the  second  central  frequency  moment  of 
the  power  spectral  density  of  the  TBL  fluctuations,  as  such  is  is  a mea- 
sure of  the  equivalent  angular-frequency  bandwidth  squared.  Furthermore, 


lim  df(u\ 


V ^^^2 

c 


27lf 


c J sin 


[?i] 


r 

L^c 


■J 


-A8- 


z. 


and 


U si 
c 


n ^ 

L _c  J cos  r 2TTfp^^”l 

Lt^J' 

c 


lim  d^ffu 


^ r T 


U 3u-^ 
c 2 


= - |(2TTfo)^  + (TTW)^  - 2^^y 

.[el 


u 

c 


cos 

L-57-J 


- 4u^fQW  cos 


cos  j 

1 1 

C =1 
n s: 

,T'l 

1 1 

FttWC" 

U 

c 

_ 

■2TTfQC 
U ^c 


] 


47Tf -U 
0 c 


sin  1 

ttW^ 
^c  . 

' ttWi 

U 

c 

sin 


2TTfoS 


cos 


TTWg 


1 

r 

c J cos 

r2TTf.n 

ttWi 

r 

U 

U 

2, 

u 

L C J 

c 

. 

2tiWU 
c 

c 


These  equations  are  combined  in  all  their  glory  to  provide 
a general  expression  for  the  variance  of  the  estimator  error  in  §3.0, 
As  a checkpoint,  the  case  for  = 0 (i.e.,  infinite  SNR)  is 

examined.  Then, 
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The  interpretation  of  the  latter  form  is  that  the  rms  error  varies 
inversely  as  the  square  root  of  the  number  of  independent  Scimples  used 
in  assembling  an  estimate  of  U^.  It  decreases  as  the  equivalent 
angular-frequency  bandwidth 

^2TTfQy  + (1/3)  (TiW) 

increases.  It  further  decreases  initially  with  an  increase  in  the 
hydrophone  separation , ^ , but  eventually  reaches  a minimum » its  decrease 
being  overta)cen  by  the  factor 


1+  exp 


r ^^0 ' 

L u I 


In  §3.0  it  is  shown  that  the  minimum  occurs  when 


5 = 1.584  U . 

^ c ^opt. 


The  rms  error  also  is  proportional  to  the  actual  value  of  the 
convection  velocity,  U^,  squared.  However,  one  of  the  U^'s  in  the 
squared  quantity  matches  dimensionally  with  the  product  of  C and  the 
equivalent  bandwidth.  When  the  optimum  separation,  is  used, 

one  of  these  U^'s  cancels,  and  the  rms  error  then  is  directly 
proportional  to  U^. 

If  the  band  in  which  the  TBL  pressure  fluctuations  are 
observed  is  too  narrow,  competing  secondary  pea)cs  evince  themselves 
in  the  space-time  correlation  function.  These  secondary  peaks  cause 
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In 


ambiguities  in  locating  the  peak  of  the  correlation,  which  the 

MLE  indicates  to  be  the  optimum  way  in  which  to  estimate  U . 

c 

order  to  avoid  this  undesirable  effect,  the  observation  frequency  band 
should  be  sufficiently  wide.  In  §3.0  a rationale  is  constructed 
for  choosing  a value  of  bandwidth  ratio  £ q/^  equal  to  3/4. 


-All- 


A. 2 Deviation  of  the  Bias  Error 


From  the  development  in  §3.0,  the  bias  error  was  given  by 


Using  the  equations  developed  in  Appendix  A.l, 


This  may  be  written 


As  m^y  be  seen,  the  bias  error  vanishes  as  the  SNR  becomes 
large,  and  the  estimator  becomes  unbiased.  For  the  values  of  y and  N 
derived  from  the  equations  in  §3.0,  the  bias  error  fraction  b ^^c)/^c 
becomes  negligibly  small.  In  fact,  once  a value  of  N is  chosen  to 
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make  the  variance  of  the  estimator  (i.e.,  fluctuation)  error 
acceptably  low,  this  value  usually  is  more  than  adequate  for  SNR's 
in  the  order  of  100  or  more  to  guarantee  a bias  error  negligible 
in  comparison  with  the  fluctuation  error. 
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A. 3 Cramer-Rao  Bound  on  Convection  Velocity  MLE 

The  Cramer-Rao  (CR)  bound  provides  a lower  bound  on  the 
minimum  variance  of  the  estimator  error.  Although  not  always 
achieved,  it  serves  as  a measure  of  the  ultimate  in  accuracy  that 
may  be  achieved  under  an  MLE.  When  the  achievable  variance  of  the 
estimator  error  is  compared  with  the  CR  bound,  a measure  of  the 
efficiency  of  the  estimator  is  derived.  In  some  cases,  the  CR  is 
achieved  as  the  number  of  samples,  N,  tends  to  infinity,  i.n  which  case 
the  bound  is  said  to  be  asymptotically  efficient. 


In  the  notation  of  §3.0,  the  CR  bound  may  be  written  as 


or 


Var 


where  E{  • } connotes  the  expected  value  of  the  quantity  within  the 
braces,  and  is  the  bias-error  given  by  ^c  ‘ 

an  unbiased  estimator,  is  zero. 


Returning 
logarithum  of  the 


to  §3.0,  the  partial  derivative  of  the  natural 
likelihood  function  reduces  to 
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As  might  be  anticipated,  the  expected  value  of  the  derivative 
vanishes.  However,  this  is  not  so  for  the  square  of  each  side  of  the  | 
equation.  Indeed,  after  some  manipulation,  } 


I 


-A16- 


becomes  (using  the  rule  for  decomposing  the  expected  value  of  the 
fourfold  product  of  Gaussian  processe*  into  the  sum  of  cyclic 
products  of  pair-wise  products;  as  exoloited  in  Appendix  A. 1) 


E 


loSe 


When  Y "*■0  (i.e.  , asymptotically  high  SNR)  , Pj^-*-Pp.  The 
CR  bound  is  derived  without  regard  for  the  data  manipulations 
indicated  by  a MLE  solution.  In  fact,  the  CR  bound  essentially 
provides  that, under  certain  conditions  of  regularity,  a greatest 
lower  bound  exists.  Other  more  stringent  conditions  such  as 
sufficiency,  and  a tighter  condition  mentioned  in  §3.0»  must  be 
observed  if  the  equality  of  the  bound  is  to  be  satisfied.  Thus,  in  view 

/ p \ (A) 

of  the  fact  that  the  observables  X and  X are  taken  initially  with- 
out knowledge  that  an  MLE  of  would  result  from  deriving  an  estimate 

of  the  time  delay  t=^  which  maximizes  the  TBL  space-time  correlation 

U 

c 

function,  the  function  pp(5,0)  is  the  one  appropriate  for  deriving 
the  CR  bound. 


Then,  for  y~'-*-0,  using  p (C,0)  and  ^Pp  , and  performing 

9U 

c 

the  indicated  manipulations,  a family  of  CR  bounds  may  be  derived 

for  each  value  of  the  variable  C-  In  fact,  it  is  readily  seen  that 

the  CR  bound  thus  derived  is  not  uniformly  the  greatest  lower  bound 

over  the  range  of  ^elO,  ®).  A minimum  value  of  the  CR  bound  is  found 

to  correspond  to  ^^0  = 0.1;  which  differs  markedly  from  the  value 

U 
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^opt^O  = 1.584  appropriate  to  the  minimum  variance  of  the  estimator 
U 

c 

error  derived  by  the  MLE  procedure.  Moreover,  the  value  of  the  CR 

bound  derived  for  ^^0  = 0.1  is  found  to  be 
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var  f U -U  » 0.  7602  U ^ 

\ c cJ  c_ 
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as  opposed  to 

var  ( U -U  ) = 0.0896  U ^ 

\ c c/  c_ 
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corresponding  to  the  MLE  for  a choice  of  ^opt^O  = 1.584.  Thus,  the 

U 

c 

MLE  procedure  has  resulted  in  a relatively  unbiased  estimate  of 
"abnormally  high"  precision;  exceeding  that  indicated  by  the  CR 
bound! 


This  anomaly  has  been  carefully  checked,  and  a lot  of  time  and 
effort  was  expended  in  chasing  an  illusory  factor  of  2n*,  but  the 
discrepancy  appears  to  be  real.  Its  roots  most  likely  are  involved 


* A factor  of  2tt  , as  either  a divisor  of  the  CR  bound  or  a multiplier 
of  the  MLE  variance,  would  make  the  MLE  of  convection  velocity  an 
estimator  of  efficiency  2;  i.e.,  the  variance  of  the  MLE  would  be 
twice  the  CR  bound.  As  observed  in  §3.0,  this  MLE  would  not  lead  to  a 
sufficient  estimator,  and  an  efficiency  of  2 might  be  expected. 
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in  the  non-regularity  of  the  PDF  of  the  observables  around  the 
value  of  ^ = 0 or  = «>.  Indeed  the  Bakewell  (loc.  cit.)  space- 
time correlation  function  has  a discontinuity  in  the  immediate 
neighborhood  of  these  values.  In  fact,  the  function 


p-0.  7 

^^0 

L 

u 

c J 

entering  into  Bakewell 's  space-time  correlation  function  corresponds 

to  the  autocorrelation  function  for  a Ist-order  Markoffian  (i.e., 

autoregressive)  process  when  i • Swerling  (loc.  cit.),  and  his 

U 

c 

doctoral  students,  investigated  other  tighter  bounds  and  estimation 
procedures  for  non-regular  estimation  cases;  including  cases  concerned 
in  estimating  correlation  coefficients  and  functions  associ- 
ated with  seemingly  innocuous  problems  involving  Ist-order 
Markoffian  processes.  This  work  involved  the  Barankin  approach;  see 
Swerling  loc.  cit.  and  Glove  loc.  cit.  In  addition,  Cramer  (loc.  cit, 
page  485)  alludes  to  similar  anomalies. 

Other  estimation  procedures  and  bounds  have  been  studied, 
such  as  the  c-a  estimator  by  Mohajeri  [6]  and  the  Ziv-Zakai  bound 
discussed  in  a paper  by  Seidman  17].  Each  claims  a different 
advantage  or  tighter  bound  under  various  assumptions.  Nevertheless, 
the  MLE  derived  in  §3.0  and  Appendix  A. 1 does  have  all  the  attributes 
intuitively  expected  of  an  estimator  of  convection  velocity.  Hence, 
for  the  purposes  of  deriving  design  equations  and  guidelines,  the 
MLE  is  offered  as  a reasonable  estimation  procedure.  It  may  not  be 
a locally  "best"  unbiased  estimator  in  the  sense  of  the  recursion 
techniques  provided  by  the  Barankin  approach;  but  sufficient  to  the 
day  is  the  rigor  thereof. 
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maximum- likelihood  estimate  of  convection  velocity  (which  is  re- 
lated to  tow  speed)  is  achieved  by  performing  a cross  correlation 
between  the  outputs  of  two  streamwise-separated  flush-mounted 
surface  probe  hydrophones.  The  estimate,  then,  is  derived  by 
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seeking  the  time  delay  which  maximizes  the  cross  correlation. 

Conditions  are  examined  consistent  with  achieving  accurate  speed 
determination  to  within  a fraction  of  a knot.  Fortuitously, 
the  method  derived  has  the  property  that,  as  speed  is  increased, 
both  the  level  and  bandwidth  of  the  pressure  signals  increase, 
thus  providing  increased  signal-to-noise  ratio  and  more  accurate 
time-delay  information. 

It  is  suggested  that  a confidence- level  measurement,  comparing 
the  space-time  correlation  function  derived  from  the  data  with 
a theoretical  model,  could  be  used  to  indicate  flow  anomalies. 
These  flow  anomalies  might  be  attributable  to  disturbances  such 
as  flow  separation.  In  such  case,  the  device  could  b«  used 
also  to  sense  flow  stability. 

Two  implementations  are  discussed:  one  using  a DElay  Line  Time 
Compressor  (DELTIC)  and  serial  processing;  the  other  using 
Shift  Registers  CSR)  and  parallel  processing.  A trade-off  of 
serial  f aster-than-real-time  DELTIC  processing  speed  require- 
ments against  the  high  degree  of  parallel  processing  needed  with 
the  SR  approach  is  discussed. 

A recommendation  is  made  to  design,  build,  and  test  such  a 
speed  estimation/stability  sensor  on  standard  towed  arrays.  It 
further  is  suggested  that  empirical  results  and  experience  thus 
derived  could  be  extrapolated,  along  with  a study  of  available 
other  array  surface  pressure  coherency  data,  to  determine  the 
feasibility  of  building  such  a sensor  for  non-standard  towed 
array  application. 
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